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Foreword 


The  Naval  Ocean  Research  and  Development  Activity  has  been  investigating 
a  possible  application  of  the  airborne  electromagnetic  method  to  bathymetric 
charting  in  a  shallow  ocean.  There  is  a  strong  Navy  requirement  for  a  rapid, 
airborne,  shallow-ocean  bathymetric  measurement  method  that  will  supple¬ 
ment  or  even  replace  the  traditional  shipbome  acoustic  sounding  methods 
that  are  timeoonsuming  and  often  not  suited  to  shallow  costal  areas.  Periodical 
and  repetitive  bathymetric  mapping  of  heavily  trafficked  shallow-ocean  regions 
is  necessary  for  meeting  specific  fleet  requirements  and  for  monitoring  bot¬ 
tom  sediment  movements,  ship  lane  maintenance,  and  a  variety  of  geotechnical 
operations,  as  well  as  for  routine  charting. 


R.  R  Onorati,  Captain,  USN 
Commanding  Officer,  NORDA 


Executive  summary 


An  experimental  airborne  electromagnetic  (AEM)  survey  was  carried  out 
in  the  Cape  Cod  Bay  area  to  investigate  the  potential  of  extracting  bathymetric 
information  for  a  shallow  ocean.  A  commercially  available  Dighem  III  AEM 
system  was  used  for  the  survey  without  any  significant  modification.  The 
helicopter-borne  system  operated  at  385  Hz  and  7200  Hz,  both  in  a  horizontal 
coplanar  configuration.  A  concurrent  ground  truth  survey  included  exten¬ 
sive  acoustic  soundings,  as  well  as  spot  water  conductivity  measurements. 

Because  of  a  lack  of  knowledge  about  the  absolute  system  calibration  figures, 
an  acoustic-sounding  calibration  was  made  for  each  flight  line  using  a  small 
portion  of  AEM  data  to  derive  the  zero-level  signal,  amplitude,  and  phase 
calibration  factors  for  each  coil  pair.  The  interpreted  bathymetric  profiles 
show  excellent  agreement  with  corresponding  acoustic  depth  profiles  up  to 
one  (possibly  more)  skin  depth  of  the  source  frequency.  It  is  envisioned  that 
with  further  improvements  in  hardware  and  software,  the  bathymetric  resolu¬ 
tion  may  extend  beyond  the  skin  depth.  AEM  data  can  also  produce  (as  by¬ 
products)  conductivity  profiles  of  both  seawater  and  bottom  sediments  that 
may  find  potential  applications  in  mine  warfare  and  offshore  geotechnical 
engineering  works. 
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Airborne  electromagnetic  bathymetry 


Introduction 

The  Naval  Ocean  Research  and  Development  Activity 
(NORDA)  has  been  investigating  a  possible  application 
of  the  airborne  electromagnetic  (AEM)  method  to 
bathymetric  charting  in  a  shallow  ocean.  There  is  a  strong 
Navy  requirement  for  a  rapid  airborne  and  cost-effective 
shallow-ocean  bathymetric  method  capable  of  supplement¬ 
ing  or  even  replacing  the  traditional  shiph  *  ne  acoustic 
sounding  methods,  which  are  time-consuming  and  often 
not  suited  to  shallow  coastal  areas.  Periodical  and  repetitive 
bathymetric  mapping  of  heavily  trafficked  shallow-ocean 
regions  is  necessary  for  monitoring  bottom  sediment 
movements,  ship  lane  maintenance,  and  a  variety  of 
geotechnical  operations,  as  well  as  for  routine  charting. 


Test  survey  in  Cape  Cod  Bay 

The  test  survey  area  and  the  AEM  flight  lines  are  shown 
in  Figure  1.  All  flights  and  ground  truth  surveys  were 
performed  during  a  3-day  period  in  June  1984.  The  AEM 
system  used  was  a  commercially  available  Dighem  III, 
described  in  detail  by  Fraser  (1978, 1979,  and  1981).  The 
system  was  equipped  with  two  horizontal  coplanar  coil 
pairs  operating  at  385  Hz  and  7200  Hz.  Both  pairs  had 
an  8-m  coil  separation  (an  additional  coaxial  coil  pair 
operating  at  900  Hz  was  deactivated  due  to  an  electronic 
malfunction). 

The  sensor  platform,  or  bird,  towed  by  a  Sikorsky  S58T 
twin-engine  helicopter  using  a  30-m  cable,  maintained  an 
average  altitude  range  between  40  m  and  50  m  above  the 
sea  surface.  The  aircraft  altitude  was  measured  by  a  radar 
altimeter  (Sperry  Model  220  mounted  on  the  aircraft)  that 
had  a  manufacturer -specified  accuracy  of  5  % .  A  total  of 
about  200  line-km  AEM  data  consisting  of  13  segment 
profiles  was  obtained  in  three  sorties  in  less  than  7  hours. 
The  sampling  rate  was  1  sec,  corresponding  to  about  50 
m  along  the  ground  track  (about  a  3  km/min  ground 
speed).  The  maximum  water  depth  in  the  survey  area  is 
about  40  m,  according  to  the  bathymetric  chart  (NOAA 
Chart  1 3,246). 

The  flight  plan  included  data  collection  before  and  after 
each  profile  at  an  altitude  of  about  270  m  to  calibrate  the 


zero-level  signal  of  the  receiver  coils.  In  addition,  three 
short  calibration  flights  were  made  in  a  location  about 
15  km  east  of  the  Cape,  where  the  bathymetric  chart  in¬ 
dicated  water  depth  in  excess  of  60  m.  These  data  were 
intended  to  be  used  for  determining  the  absolute  calibra¬ 
tion  constants  for  amplitude  and  phase  of  each  coil  pair 
on  an  assumption  that  the  water  body  below  may  be  con¬ 
sidered  a  uniform  conductive  half-space.  It  turned  out, 
however,  that  this  calibration  method  is  not  accurate 
enough  for  the  bathymetric  processing.  As  discussed  in 
“Data  Calibration,”  both  zero-level  signal  and 
amplitude/phase  calibration  constants  are  derived  from  a 
small  portion  of  each  actual  flight  line  data. 

Figure  2  shows  a  raw  AEM  data  profile  accompanied 
by  a  corresponding  radar  altimeter  profile  along  Line  5021 
(see  Fig.  1  for  location).  Clearly,  the  AEM  data  are  over¬ 
whelmingly  correlatable  with  variations  in  altitude.  A  very 
crude  indication  of  water  depth  may  be  observed  from 
the  ratio  of  the  quadrature  component  to  the  inphase  com¬ 
ponent  of  the  385-Hz  data:  the  ratio  increases  with  a 
decreasing  water  depth.  Unfortunately,  this  relationship 
is  highly  nonlinear.  Even  though  the  aircraft  altitude  is 
maintained  mostly  within  a  10-m  range  (between  40  m 
and  50  m),  the  corresponding  variations  of  the  AEM 
responses  amount  to  more  than  500  parts  per  million 
(ppm).  Owing  to  the  high  water  conductivity,  errors  in¬ 
duced  by  inaccurate  altimetry  pose  a  critical  problem.  At 
a  45-m  bird  altitude,  a  1  %  altitude  change  at  a  given  water 
depth  of  10  m  generates  amplitude  differences  of  22  ppm 
at  385  Hz  and  33  ppm  at  7200  Hz.  It  can  also  be  shown 
that,  for  a  1-m  depth  change  at  the  same  water  depth  of 
10  m,  the  predicted  amplitude  differences  amount  to 
only  10  ppm  at  385  Hz  and  0  ppm  at  7200  Hz. 

Since  the  employed  radar  altimeter  has  a  specified  ac¬ 
curacy  of  5%,  it  soon  became  evident  that  the  radar 
altitude  cannot  be  trusted  for  the  bathymetric  processing. 
Instead,  a  new  algorithm  was  developed  to  use  the  7200 
Hz  response  to  derive  the  electromagnetic  altitude  dur¬ 
ing  the  inversion  process.  The  new  altitudes  thus  derived 
show  fairly  random  zero-biased  differences  (with  respect 
to  the  radar  altitudes)  whose  rms  amounts  to  about  2-3%. 

Navigation  was  originally  planned  to  employ  a  Del 
Norte  navigation  system  supported  by  three  ground 
transponders.  Excessive  distances  caused  p<x>r  reception: 
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AEMB  LINE:  3021 


Figure  2.  Raw  AEM  and  radar  altitude  data  along  Flight  Line  5021. 


therefore,  a  Loran-C  system  was  installed  on  site  with  a 
makeshift  arrangement  of  a  printer  that  produced  coor¬ 
dinates  at  a  5-sec  interval.  These  were  later  interpolated 
to  produce  1-sec  interval  coordinate  data  corresponding 
to  the  AEM  data  rate. 

A  ground  truth  bathymetric  survey  concurrent  with 
the  AEM  flights  was  carried  out  using  an  acoustic  depth 
sounder.  A  total  of  about  120  line-km  depth  profiles  was 
obtained,  which  covered  about  60%  of  the  AEM  flight 
area.  Unfortunately,  due  to  many  practical  reasons,  the 
flight  lines  and  the  ship  track  did  not  coincide  and  were 
often  more  than  500  m  apart.  Therefore,  the  best  available 
ground  truth  still  reflects  another  interpolated  approxima¬ 
tion  (unless  the  bottom  topography  fluctuates  rapidly,  the 
ground  truth  is  considered  to  be  accurate  within  1  to  2  m). 

Spot  measurements  of  water  conductivity  were  made 
at  eight  different  locations  along  the  ship  track  at  a  3-m 
depth.  They  ranged  between  4.0  mho/m  and  4.12  mho/m. 
While  these  values  may  be  fairly  representative  for  deep 
water,  there  are  considerable  uncertainties  over  very 
shallow  water  ( <  3  m)  where  water  temperature  may  rise 
significantly  during  the  day  (particularly  during  sunny  days 
in  June,  as  in  this  case).  A  mere  4°C  difference  in  the 
water  temperature  at  a  given  salinity  can  result  in  as  much 
as  a  10%  change  in  water  conductivity.  Unfortunately, 
no  ground  truth  measurements  were  made  during  the 
survey  to  confirm  this  possibility. 


Interpretation 

The  high  conductivity  of  seawater  (between  3  and  5 
mho/m,  depending  on  salinity  and  temperature  with  no 
fresh  water  inlets)  severely  restricts  the  ability  of  EM  waves 
to  penetrate  the  water.  Bathymetric  range  and  resolution 
are,  therefore,  primarily  governed  by  the  source  frequen¬ 
cy.  Figure  3  shows  the  skin  depths  in  a  frequency  range 
between  40  Hz  and  40  kHz  for  assumed  water  conduc¬ 
tivities  of  2,  3,  4,  and  5  mho/m. 

For  the  employed  frequencies  of  385  Hz  and  7200  Hz 
for  seawater  with  a  conductivity  of  4  mho/m,  we  may, 
therefore,  expect  skin  depth  of  12.8  m  and  3.0  m,  respec¬ 
tively.  From  Figure  3  the  source  frequency  obviously 
should  be  less  than  100  Hz  to  achieve  a  depth  range  of 
50  m  or  more. 

Fundamental  equations  for  the  magnetic  field  generated 
by  a  vertical  magnetic  dipole  located  at  or  above  the  sur¬ 
face  of  a  layered  earth  are  given  by  Kozulin  (1963)  and 
Frischkneckt  (1967).  The  mutual  coupling  ratio  for  a 
horizontal  coplanar  configuration,  used  for  the  present 
frequency-domain  AEM  system,  is  defined  as  the  ratio 
of  the  total  magnetic  field  (Hz)  to  the  primary  field  (f/f): 

OD 

Hz/Hf  =  1  -  a3  j  \2  R  (X,  d,  at,  a ?  h,f)  J0  (Kj)  cfK . 
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Figure  3-  Skin  depths  as  a  function  of  frequency  for  seawater  having  a  conductivity  ranging  from  2  mho/m  to  5  mho/m. 


The  kernel  function  R  corresponding  to  a  two-layer  earth, 
of  which  geometry  is  shown  on  Figure  4,  can  be  expressed 
as 

is  .  is  --2VV 
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R  =  - ,  (2) 
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where  VQ  ]  =  (V0  -  V,)  /  (Vn  +  V,J, 

VI  2  =  (V,  -  VJ  /  (V,  +  Vf), 

K  =  X _ 

Vt  =  V\2  +  i  2ir  p.  a,  f ,  and 

V2  =  VX2  +  i  2r  ft  a2f  ■ 

The  mathematical  notations  in  the  above  expressions  are: 
/:  transmitter  frequency  (Hz), 
h\  bird  altitude, 
a:  coil  separation, 
d\  water  depth, 

<7;:  water  conductivity, 
a2:  sediment  conductivity, 
ft:  magnetic  permeability  of  free  space, 

X:  variable  of  integration,  and 
J0:  the  zeroth  order  Bessel  function. 
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Figure  4.  The  A  EM  bathymetric  model  used  for  inversion 
The  unknowns  are  water  depth,  conductivities  of  seawater 
and  bottom  sediment,  and  the  bird  altitude. 
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The  integral  in  Equation  (1)  can  be  evaluated  by  the  linear 
digital  filter  method  (Koefoed  et  al.,  1972).  We  used  the 
filter  coefficients  published  by  Anderson  (1979b)  for  the 
Hankel  transform  integral. 

The  first  term  in  Equation  (1)  representing  the  primary 
field  is  customarily  bucked  out  during  measurements,  and 
only  the  second  term  representing  the  ocean  response  is 
recorded  in  a  ppm  unit.  Figure  5  shows  computed  inphase 
and  quadrature  frequency  responses  in  a  range  of  40  Hz 
to  40  kHz  for  various  water  depths  up  to  50  m  for  a  50-m 
bird  altitude,  an  8-m  coil  separation,  and  conductivities 
of  water  and  sediment  4  mho/m  and  1  mho/m,  respectively. 

Compared  with  normal  land  survey,  the  ppm  responses 
over  an  ocean  are  extremely  high  due  to  the  highly  con¬ 
ductive  seawater.  It  is  also  noted  that  the  responses  are 
critically  sensitive  to  the  bird  altitude.  The  inphase  response 
increases  with  the  source  frequency  while  the  quadrature 
response  peaks  approximately  at  a  frequency  at  which  the 
water  depth  equals  the  corresponding  skin  depth.  This 
is  previously  explained  theoretically  by  Won  (1980). 

Since  we  are  concerned  with  the  differential  changes 
in  the  response  for  varying  bathymetric  depth,  we  pre¬ 
sent  Figure  6,  which  shows  the  differences  of  the  ppm 
responses  with  respect  to  an  infinitely  deep  ocean  having 
the  same  4  mho/m  conductivity.  We  now  notice  that  the 


total  change  corresponding  to  a  water  depth  change  from 
10  m  to  50  m  amounts  to  about  100  ppm.  In  general, 
this  response  increases  proportionally  to  the  inverse  cube 
of  the  bird  altitude.  It  is  obvious,  therefore,  that  the  record¬ 
ing  device  must  have  much  larger  dynamic  ranges  than 
that  used  for  land  survey. 

Various  inversion  techniques  are  presently  available  to 
solve  for  the  unknown  parameters  in  Equation  (1).  These 
include  predominantly  several  variations  of  the  least-squares 
method,  including  the  Marquardt  algorithm  (Marquardt, 
1963;  Anderson,  197%)  and,  occasionally,  applications  of 
the  generalized  inverse  theory  (Backus  and  Gilbert,  1967; 
Fullagar  and  Oldenburg,  1984;  Son,  1985). 

The  Cape  Cod  test  data  were  initially  interpreted  and 
reported  by  Fraser  (1985)  using  a  least-squares  algorithm 
by  Anderson  (1979a).  Subsequently,  the  data  were 
reprocessed  at  NORDA  using  a  different  Marquardt  least- 
squares  algorithm,  notably  Subroutine  ZXSSQ  in  the  IMSL 
package.  The  inverted  bathymetry  in  both  cases  agreed 
approximately  in  trends  with  known  bathymetry  but 
showed  a  considerable  static  bias  that  often  exceeded  5-10 
m.  Further  careful  inspection  of  the  least-squares  inver¬ 
sion  results  leads  us  to  the  following  conclusions: 

•  Computer  inversion  time  is  unacceptably  long:  one- 
point  inversion  of  the  two-frequency  data  consumes 
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Figure  5  AEM  bathymetric  responses  expressed  in  ppm  for  water  depths  of  10  m,  20  m,  JO  m,  40  m  and  50  m.  The  bird 
altitude  is  assumed  to  be  50  m,  the  coil  separation  8  m,  water  conductivity  4  mho/m,  and  sediment  conductivity  I  mho/m. 


from  5  sec  to  1  min  on  a  VAX  11/780  computer,  even 
when  the  water  depth  is  the  only  sought  parameter, 
while  all  other  parameters  are  prescribed  and  fixed. 

•AEM  response  is  too  sensitive  to  the  bird  altitude 
to  accept  the  specified  5%  accuracy  of  the  radar 
altimeter  used  for  the  survey. 

•  Both  the  water  and  the  sediment  conductivities  must 
be  allowed  to  float  albeit  in  constrained  ranges. 

The  AEM  bathymetric  profiles  reported  here  are  derived 
from  yet  another  method:  analytic  solutions  of 
simultaneous  nonlinear  equations.  At  each  data  location, 
we  have  four  measured  quantities;  i.e.,  inphase  and 
quadrature  components  at  two  frequencies.  From  this  data 
set  we  derive  exact  solutions  of  four  parameters:  water 
depth,  water  conductivity,  sediment  conductivity,  and  elec¬ 
tromagnetic  altitude.  When  unconstrained,  the  solutions 
are  exact  (since  the  number  of  knowns  and  unknowns  is 
the  same),  resulting  in  zero  residuals  regardless  of  data 
error.  However,  severe  data  error  may  produce  physical¬ 
ly  unacceptable  solutions  (e.g.,  negative  depth  or  conduc¬ 
tivities).  While  least-squares  methods  (in  which  the  number 
of  knowns  is  usually  much  more  than  that  of  unknowns) 
may  produce  a  stable  solution  set  (even  though  its  rms 
error  may  be  high)  from  a  noisy  data  set,  the  present 
analytic  approach  is  understandably  sensitive  to  data 


error.  Under  this  circumstance,  a  low-pass  filtering  of  the 
inverted  profile  is  justifiable  to  countermeasure  the  ran¬ 
dom  data  errors. 

An  inversion  algorithm  using  a  modified  Newton- 
Raphson  method  is  then  applied  to  the  data.  Initially,  we 
derive  the  sensor  altitude  and  water  conductivity  from 
the  7200-Hz  data  and,  subsequently,  water  depth  and  bot¬ 
tom  conductivity  from  the  385-Hz  data.  Inversion  time 
for  deriving  all  four  parameters  amounts  to  0.5  to  2  sec 
on  a  VAX  11/780  computer.  The  analytic  method,  as 
in  the  least-squares  method,  also  requires  initial  guesses 
and,  to  ensure  physically  acceptable  solutions,  reasonable 
solution  constraints.  The  constraints  used  for  the  final  proc¬ 
essing  of  the  Cape  Cod  data  follow. 

•  Water  conductivity  (a,):  3-5  mho/m 

•  Sediment  conductivity  ( o2 ):  0.01-2  mho/m 

•  Water  depth  (d)\  0-50  m 

•  Altitude  ( h ):  positive 

Spot  measurements  of  water  conductivity  at  a  3-m  depth 
at  eight  locations  ranged  from  4.0  to  4.12  mho/m.  No 
bottom  sediment  conductivity  data  are  available.  However, 
an  extensive  in  situ  study  by  Hulbert  et  al.  (1982)  off  the 
Florida  coast  shows  a  common  range  of  0.4  mho/m  to 
1.4  mho/m  within  the  first  5-m  depth,  decreasing  only 
slightly  with  increasing  depth  of  burial. 


Figure  6  AEM  bathymetric  response  expressed  in  ppm  for  water  depths  10  m,  20  m,  30  m,  40  m  and  50  m  relative  to  fl 

an  infinitely  deep  water.  The  bird  altitude  is  assumed  to  be  50  m,  coil  separation  8  m,  water  conductivity  4  mho/m.  and  j 
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The  inversion  process  is  initiated  as  follows:  For  the 
very  first  point,  we  prescribed  starting  values  of  ar  = 
4  mho/m,  a,  =  1  mho/m,  d  as  read  from  the 
hydrographic  chart,  and  h  as  indicated  by  the  radar 
altimeter.  Once  the  process  starts,  the  solution  set  at  the 
present  location  is  prescribed  as  the  initial  parameters  for 
the  next  location.  Thus,  after  the  first  data  point  of  a  pro¬ 
file,  the  interpretation  becomes  completely  autonomous. 

We  present  only  the  bathymetric  results.  Presentation 
of  other  parameters  will  be  dealt  with  in  a  separate  report. 
It  is  noted,  however,  that  (1)  the  derived  electromagnetic 
altitude  is  well  within  ±  1  m  of  the  radar  altitude  (less 
than  the  manufacturer-specified  5%  error),  (2)  water  con¬ 
ductivity  is  mainly  4  ±  0.2  mho/m  except  for  very  shallow- 
water  regions,  and  (3)  bottom  sediment  conductivity  ranges 
between  0.5  mho/m  and  1.5  mho/m  in  most  profiles. 

Data  calibration 

From  the  beginning  it  was  realized  that  using  an  ex¬ 
isting  commercial  frequency-domain  A  EM  system  designed 
for  over-the-land  survey  posed  a  serious  problem  in 
establishing  the  zero-level  signal,  and  gain  and  phase 
calibration  factors.  Exploration  geophysicists  usually  pay 
more  attention  to  relative  anomalies  than  to  their  absolute 
values.  In  the  AEM  bathymetric  survey,  we  face  the 
challenge  of  determining  the  absolute  values. 

An  initial  attempt  to  use  a  set  of  uniform  calibrations  (de¬ 
rived  from  the  deep  ocean  data)  for  the  entire  survey  data 
produced  unsatisfactory  results:  while  the  AEM  bathyme¬ 
try  approximately  followed  known  depth  profiles,  it  mani¬ 
fested  significant  static  bias  often  amounting  to  5-10  m. 

For  the  Cape  Cod  test  data,  we  experimented  with  three 
different  calibration  techniques,  viz  ,  zero-level  calibra¬ 
tion  only,  amplitude/phase  calibration,  and  zero-level/ 
amplitude/phase  calibration.  Of  these,  the  last  approach 
turned  out  to  be  superior  to  others  in  inversion  results 
and  was  adopted  for  the  final  processing. 


Calibration  factors  for  selected  profiles  are  listed  in  Table 
1.  We  note  that  the  amplitude  correction  factors  range 
between  0  and  7%,  the  phase  correction  between  2°  and 
5°,  and  the  zero-level  correction  between  -  28  and  -  2 
ppm.  Such  insignificant  corrections  may  be  quite  ignorable 
for  many  routine  mining  exploration  problems  where  only 
relative  anomalies  are  sought. 

The  calibration  constants  are  derived  as  follows:  to  stand¬ 
ardize  the  process,  we  choose  an  arbitrary  50-data  point 
profile  segment  (about  2.4  km  long  for  the  present  data) 
over  a  relatively  flat  and  deep  ocean  where  an  average 
water  depth  is  known.  For  each  data  point,  we  then 
prescribe  Oj  =  4  mho/m,  o2  =  \  mho/m,  and  h  =  the 
radar  altitude.  The  last  is  the  only  parameter  that  varies 
for  each  data  point.  Using  these  prescribed  parameters, 
we  subject  the  entire  segment  data  set  to  the  Marquardt 
least-squares  inversion  to  derive  the  best  fit  amplitude, 
phase,  and  zero-level  values.  Because  of  the  large  data  set, 
the  inversion  produces  very  stable  calibration  factors.  These 
calibration  factors  are  applied  to  each  raw  data  profile 
before  the  final  inversion  process.  It  should  be  pointed 
out  that  this  known-depth-point  calibration  method  also 
compensates  for  the  tidal  fluctuation,  which  amounted  to 
a  maximum  height  of  2.8  m  during  the  7-hour  flight 
period. 

Obviously  such  a  hindsight  calibration  technique  is  unac¬ 
ceptable:  a  future  production  AEM  bathymetry  system 
must  contain  an  automatic  electronic  calibration  capabili¬ 
ty.  It  should  be  noted,  however,  that  the  data  used  for 
the  calibration  comprise  only  a  fractional  segment  (about 
10%)  of  a  given  profile;  thus,  most  of  the  profile  is  not 
directly  influenced  by  the  scheme. 

Results 

Figure  7  shows  the  interpreted  AEM  bathymetry  for 
Line  5021  (see  Fig.  1  for  location).  The  solid  line  represents 


Table  1 .  Cape  Cod  AEM  data  calibration  constants. 


Line 

385  Hz 

Ampl 

Phase 

(deg) 

Zero  level  (ppm) 

Inphase  Quad 

1.063 

2.56 

-17.1 

1.028 

3  68 

-  2  4 

5041 

1.072 

5.20 

9 

-27  5 

5051 

1  025 

3  49 

-  5.0 

1021 

1.040 

-0.61 

■  V 

38  6 

5082 

0.996 

2  54 

-  4.3 

7200  Hz 


Ampl 


Phase 

(deg) 


Zero  level  (ppm) 
Inphase  Quad 


1.089 
1  056 
1  099 
1.054 
1  078 


2  91 
4  30 
4  90 
4  30 
2  77 
2.39 


0 

0 

0 

0 

0 

0 


-12  7 
-13  3 
-27  7 
-17  3 
-14  0 
-  33 


1  036 
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Figure  7.  AEM  bathymetric  profile  for  Line  5021  along  with  radar  altimeter  profile.  Solid  line  represents  AEM  bathymetry; 
while  solid  circles  represent  depths  obtained  from  a  shipbome  acoustic  profiler.  Small  numerals  at  bottom  are  the  flight  fiducials. 


the  water  depth  inferred  from  the  AEM  data.  Solid  circles 
denote  the  depths  determined  from  acoustic  profiles. 
Depths  are  computed  at  approximately  50-m  intervals. 
Small  numerals  at  the  bottom  are  the  flight  line  fiducials 
representing  every  20th  data  point.  The  profile  length  is 
about  20  km. 

The  agreements  are  excellent  up  to  a  water  depth  cor¬ 
responding  to  about  one  skin  depth  (12.8  m)  of  the  385-Hz 
signal.  In  fact,  the  agreements  up  to  this  depth  are  well 
within  the  interpolation  accuracy  of  ground  truth  data. 
Below  the  skin  depth  we  notice  progressively  degrading 
resolution  resulting  in  oscillatory  bathymetric  profiles.  Such 
oscillatory  behaviors  over  deep  water 

•  are  common  to  all  AEM  bathymetric  profiles  obtained 
in  the  Cape  Cod  Bay, 

•  exhibit  more  or  less  the  same  rate  of  degradation  with 
depth, 

•  are  strongly  correlated  with  variations  in  the  aircraft 
altitude,  and 

•  appear  to  be  of  random  Gaussian  error  (not  rigorous¬ 
ly  determined). 

In  essence,  the  oscillatory  behavior  is  a  direct  result  of 
the  decreasing  signal-to-noise  ratio  with  respect  to  the 
altitude  uncertainty.  At  a  20-m  depth,  for  instance,  the 
maximum  theoretical  385-Hz  response  (against  an  infinitely 
deep  water)  is  expected  to  be  about  10  ppm  (Fig.  6),  while 
a  mere  0.2-m  error  in  altitude  will  result  in  the  same 
amount  of  difference  in  response.  Since  the  bathymetric 
errors  appear  to  be  random,  yet  strongly  correlated  with 
the  aircraft  altitude,  we  tentatively  conclude  that  the  error 


sources  are  likely  related  to  the  altimeter  resolution  and 
to  such  bird  attitude  uncertainties  as  pitching  and  yaw¬ 
ing  associated  with  the  aircraft  altitude  variations.  The 
bird  attitude  can  be  monitored  in  the  future  using  in¬ 
clinometers  whose  output  can  be  incorporated  into  the 
interpretation  (Son,  1985). 

Such  an  oscillatory  behavior  can  sometimes  be  sup¬ 
pressed  if  we  use  instead  a  least-squares  inversion  method 
when  a  sufficient  number  of  redundant  measurements  is 
available.  The  resultant  solutions  in  this  case  will  carry 
large  rms  errors,  yet  may  give  a  deceptively  smooth  solu¬ 
tion  profile  (errors  never  die;  they  simply  become  hidden 
in  the  process).  The  present  analytic  inversion  method 
produces  zero-residual  solutions  that  fit  the  observed  data 
regardless  of  the  measurement  errors.  Although  the  two 
approaches  are  equivalent  in  the  sense  of  error  budgeting, 
the  analytic  inversion  method  appears  to  be  superior  in 
field  logistics  and  in  computational  speed. 

Figure  8  shows  the  interpretive  error  profiles  for  Line 
5021.  The  top  figure  shows  the  difference  between  the 
radar  altitude  and  the  electromagnetic  altitude  derived  from 
the  7200-Hz  record.  The  differences  are  mainly  less  than 
1  m,  except  toward  the  shoreline  where  water  becomes 
very  shallow  and  where  the  interpretative  model  for  the 
AEM  data  fails.  While  the  differences  appear  to  be  ran¬ 
dom,  they  are  fairly  correlatable  with  the  altitude  changes. 
It  is  understandable  that  when  aircraft  altitude  changes, 
the  air  speed  in  general  also  changes,  thus  resulting  in 
changes  in  the  towing  angle  which,  in  turn,  produces  a 
relative  motion  between  the  aircraft  (on  which  the  altimeter 


is  mounted)  and  the  bird.  This  relative  motion  appears 
to  be  responsible  for  the  differences.  Mounting  the 
altimeter  on  the  bird  may  help  to  resolve  this  problem 
in  the  future. 

If  we  assume  that  the  oscillatory  behavior  of  the  AEM 
bathymetry  (Fig.  7)  is  of  random  nature,  we  are  justified 


to  perform  a  low-pass  filtering  of  the  interpreted 
bathymetric  profile  to  render  smooth  appearances.  To  this 
end,  we  applied  to  Line  5021  a  simple,  equal-weight, 
11-point  running  average  filter  to  produce  Figure  9.  Figures 
10-14  show  additional  AEM  bathymetric  profiles  produced 
by  the  above  described  procedure.  The  same  11 -point  filter 
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Figure  8.  Interpretative  error  profiles  for  Line  5021.  Top  graph  shows  the  differences  between  the  radar  altitude  and  the 
electromagnetic  altitude  deduced  from  the  7200  Hz  data.  The  rest  shows  differences  between  measured  EM  responses  and 
computed  EM  responses  for  385  Hz  and  7200  Hz  records. 
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Figure  9  AEM  bathymetric  profile  for  Line  5021  after  applying  an  11  point  running  average  filter.  Solid  circles  represent 
acoustic  depths. 
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AEMB  LINE:  5031 


Figure  10.  AEM  bathymetry  profile  for  Line  5031.  Solid  circles  represent  acoustic  depths. 


AEMB  LINE:  5041 


has  been  applied  to  all  profiles.  Where  the  ground  truth 
survey  was  not  performed,  we  show  water  depths  as  read 
from  the  bathymetric  chart. 

A  composite  of  seven  AEM  profiles  is  shown  in  Figure 
15.  We  notice  striking  details  of  the  sea  bottom  mor¬ 
phology  showing  subtle  trends  and  developments  of  slopes, 
trenches,  and  shoals.  The  fact  that  each  profile  is  in¬ 
dependently  derived  and  yet  shows  remarkable  correla¬ 
tions  with  neighboring  profiles  renders  further  credence 
to  the  AEM  results. 


Conclusions 

From  our  experience  through  the  Cape  Cod  AEM 
bathymetry  experiment,  we  summarize  some  of  the  er¬ 
ror  sources  that  degrade  the  bathymetric  resolution: 

•  calibration  errors:  amplitude,  phase  and  zero-level; 

•  error  in  the  interpretative  ocean  model,  particularly 
assuming  the  vertically  homogeneous  bottom  sediment 
layer; 

•altimeter  error; 
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Figure  12.  AEM  bathymetry  profile  for  Line  5051.  Solid  circles  represent  acoustic  depths. 
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Figure  13.  AEM  bathymetry  profile  for  Line  1021.  Solid  circles  represent  acoustic  depths. 


•  measurement  error  due  to  pitching  and  yawing  of  the 
bird— negligible  up  to  10°  if  the  bird  altitude  is  50 
m  or  higher; 

•  ground  truth  interpolation  error  due  to  noncoincidence 
of  tracks  by  boat  and  aircraft; 

•  electronic  measurement  noise. 

Most  of  the  above  error  sources  can  be  significantly  re¬ 
duced  through  improvements  in  equipment  and  interpreta¬ 
tion  software. 

It  is  envisioned  that  with  additional  research  and  develop¬ 
ment  efforts,  the  AEM  method  will  be  able  to  produce 


accurate  bathymetric  charts  over  a  shallow  ocean  (perhaps 
up  to  100  m  in  depth).  Compared  with  the  traditional 
acoustic  sounding  techniques,  the  AEM  method  can  pro¬ 
vide  an  order-of-magnitude  faster  survey  speed  at  a  reduced 
cost  and  thus  yield  a  synoptic  knowledge  of  ocean-bottom 
topography.  With  improved  interpretation  schemes,  even 
real-time  data  processing  appears  to  be  a  realizable  goal. 

In  addition,  the  method  has  potential  applications  to 
remote  measurements  of  electrical  conductivities  of  ocean 
water  and  bottom  sediments.  The  bottom  sediment  con¬ 
ductivity,  in  particular,  is  closely  related  to  certain 
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Figure  14 ■  AEM  bathymetry  profile  for  Line  5082.  Solid  circles  represent  water  depths  read  from  the  bathymetric  chart  (NO A  A 
chart :  13,246). 
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Figure  15.  A  composite  of  seven  AEM  bathymetry  profiles  from  the  Cape  Cod  Bay. 


mechanical  characteristics,  such  as  compaction  rate,  porosi¬ 
ty.  density,  and  (indirectly  perhaps)  sediment  types,  which 
carry  broad  geotechnical  implications  for  many  offshore 
activities. 
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